Finding the Higgs Boson through Supersymmetry 
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The study of displaced vertices containing two b-jets may provide a double discovery at the Large 
Hadron Collider (LHC): we show how it may not only reveal evidence for supersymmetry, but also 
provide a way to uncover the Higgs boson necessary in the formulation of the electroweak theory 
in a large region of the parameter space. We quantify this explicitly using the simplest minimal 
supergravity model with bilinear breaking of R-parity, which accounts for the observed pattern of 
neutrino masses and mixings seen in neutrino oscillation experiments. 
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I. INTRODUCTION 

By opening the exploration of the new territory of 
physics at the Tcrascale, the CERN Large Hadron Col- 
lider (LHC) is likely to shed light upon the main open 
puzzle in particle physics, namely the origin of mass and 
the nature of electroweak symmetry breaking. Super- 
symmetry (SUSY) provides an elegant way of justifying 
the electroweak symmetry breaking mechanism in terms 
of an elementary Higgs particle, alleviating the so called 
hierarchy problem [l| . The Higgs boson and the existence 
of supersymmetry therefore stand out as the main miss- 
ing pieces in our understanding of fundamental forces, 
and a lot of effort has been put into their direct obser- 
vation. Indeed the search for the Higgs boson and for 
supersymmetry constitute the main topic in the agenda 
of the LHC. 

In contrast, so far the only established evidence for 
physics beyond the standard model (SM) has been the 
discovery of neutrino masses and oscillations which 
has culminated decades of painstaking efforts. 

Here we stress that these two issues may be closely 
related. Indeed, low-energy supersymmetry with broken 
R-parity 0] provides a plausible mechanism for the ori- 



gin for neutrino masses and mixings. Indeed, as the bilin- 
ear model best illustrates in contrast to the simplest 
seesaw schemes 0| , these may be tested at particle accel- 
erators like the LHC 1 . 

Here we consider the simplest ansatz to introduce R- 
parity breaking in supersymmetry, characterized by an 
additional bilinear violating (BRpV) term in the super- 
potential It provides the simplest effective descrip- 
tion of a more complete picture containing additional 
neutral heavy lepton lfj superfields whose scalars drive 
the spontaneous breaking of R-parity 11 1 . 



Our focus here is on the specific case of a minimal grav- 
ity mediated supersymmetry breaking model with bilin- 
ear R parity violation: BRpV-mSUGRA model for short. 
In this model, the lightest supersymmetric particle (LSP) 
is no longer stable. Current neutrino oscillation data in- 
dicate that the strength of the BRpV term is small 0], 
hence the LSP decay length is expected to be l ong enough 



to provide a displaced vertex at the LHC [12|, |13| . For a 
low Higgs mass the dominant decay is into bb, however at 
the LHC the overwhelming QCD background makes this 
signal irrelevant when the Higgs is produced in the stan- 
dard way. In supersymmetry the Higgs can be produced 
after the decay chains of the next-to-lightest supersym- 
metric particle. In the R-parity conserving case for spe- 
cific spectrum and supersymmetric production, the addi- 
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1 Such model has no conventional neutralino dark matter, though 
other possible dark matter candidates may be envisaged such as 
the axion [f|, the majoron Q, the axino or the gravitino @. 
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tional jets and the missing energy can allow the discovery 
of the Higgs in the b channel 14[ . The same features also 
hold in our case, but in addition now the Higgs can be 
produced from the lightest ncutralino. leading to events 
with a displaced vertex with two large invariant mass 
b-jets. The signal of a neutralino into a Higgs and a 
neutrino is therefore free of SM backgrounds if the neu- 
tralino decays inside the pixel detector and well outside 
the interaction point. Here we show explicitly that this 
is the case 2 . 

In this work we analyze the potential of the LHC to 
survey the existence of the Higgs boson using a novel sig- 
nal: a b-jet pair coming from displaced vertices gener- 
ated by the lightest neutralino decays within the BRpV- 
mSUGRA model. We demonstrate that the LHC reach 
is capable to uncover a supersymmctric Higgs in a fair 
region of the M 1 / 2 <8> Mq parameter plane. 



II. MODEL DESCRIPTION 



The BRpV model is described by the superpotential 



W^BRpV = W^MSSM + £ab£iLiH^ , 



(1) 



in which the standard minimal supersymmctric model 
(MSSM) is supplemented by three extra bilinear terms 
characterized by three new parameters (ej), one for each 
fermion generation. In addition to these we must also 
include new soft supersymmetry breaking terms (Bi) in 
whose presence the bilinears become physical parameters 
that can not be rotated away fir] ]. 



Vsott = Vmssm — £ a bBieiL°;H^ 



(2) 



The new terms in the BRpV Lagrangian (e;, £?j) lead 
to the explicit violation of lepton number as well as R- 
parity. Furthermore, the sneutrino fields acquire a vac- 
uum expectation value when we minimize the scalar po- 
tential. 

In BRpV models the terms presenting explicit lepton 
number violation, as well as, the sneutrino vacuum expec- 
tation values generate mixing among neutrinos and neu- 
tralinos giving rise to one tree-level neutrino mass. The 
other two neutrino masses are generated through loop 



diagrams [9(. One can show that, indeed, the resulting 
neutrino masses and mixings provide a good description 
of all current neutrino oscillation data [2j. 

For the sake of definiteness, we assume mSUGRA as 
the model of supersymmetry breaking, implying univer- 
sality of the soft breaking terms at unification. In this 
case, our model depends upon eleven free parameters, 
namely 

M , M 1/2 , tan/3 , sign(^) , A , , and A, , (3) 

where Mi/ 2 and Mq are the common gaugino mass and 
scalar soft SUSY breaking masses at the unification scale, 
Aq is the common trilinear term, and tan/3 is the ratio 
between the Higgs field vev's. For convenience, we trade 
the soft parameters Bi by Aj = ejW<j + /iu^, where Vi 
is the vacuum expectation value of the sneutrino fields, 
since the A^'s are more directly related to the neutrino 
masses; for further details see [9(. 

The bilinear R-parity violating interaction gives rise 
to mixings between SM and SUSY particles that lead to 
decay of the LSP into SM particles. In a large fraction 
of the parameter space the lightest neutralino is the LSP 
and it can decay into leptonic final states v£ + £'~ , where 
£ = e, fi or r, as well as into scmi-lcptonic final states £q'q 
or vqq. For sufficiently heavy neutralinos these decays are 
dominated by two-body channels like vZ, £ ± 1U T and vh 
with h bein g the lig htest CP even Higgs boson; for further 

details sec El El El- in the region where the stau is 

the LSP the detached vertex signal disappears completely 
since the stau possesses a very small decay length. 

In contrast, a salient feature of our BRpV model is 
that neutralino LSPs exhibit a rather large decay length, 
ranging from a few millimeters to tenths of millimeters 
for Mx/2 varying from 200 GeV to 1 TeV. Such large 
decay lengths lead to the production of detached vertices 
at the LHC which constitute a smoking gun of this kind 
of models. 

In this work, we analyze the two-body lightest neu- 
tralino decay into the lightest Higgs boson h° as a Higgs 
discovery channel 



X? 



hv 



(4) 



2 In fact the LHCb collaboration is considering the possibility to 
search for b's originating outside the interaction point |15H . 



If the lightest neutralino lives long enough it will be de- 
tached from the primary interaction point leaving a dis- 
placed vertex as signal at the LHC. Since the Higgs bo- 
son h decays mostly into a b-quark pairs we expect a 
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Figure 1: Br(x? -> as a function of Mi/2®Mo for tan/? 
10, A = -100 GeV and fi > 0. 



displaced vertex with two b-jets as a characteristic sig- 
nature for Higgs production. 

We present, in Figure[TJ the lightest ncutralino branch- 
ing ratio to hv as a function of Mj/2 ® Mo for tan j3 = 10, 
A = —100 GeV and /i > 3 . Here we focus on the situa- 
tion where the lightest neutralino is heavier then h, so the 
neutralino Higgs decay channel opens for M1/2 > O(300) 
GeV for our choice of parameters. The maximum value 
of the branching ratio for this channel is about 22%; for 
an illustration of the full behavior of neutralino decays 
see, for example, Ref. 13, 17, [3]. This figure tells us 
that, for fixed values of M1/2, the LSP branching ratio 
into Higgs neutrino pairs initially grows with increasing 
Mo, stabilizing for Mo in excess of a few hundred GeV. 
On the other hand, the importance of this decay increases 
with M 1 j 2 for moderate and large values of Mq. 



III. SIGNAL AND BACKGROUNDS 

In order to simulate the Higgs production we calculate 
all R-parity violating branchin g ra tios and SUSY spec- 
tra using the package SPheno Hi]. We used PYTHIA 
version 6.408 [20j to generate events, using the SPheno 
output in the SLHA format [2]} . In order to have a rough 
simulation of the detector response we smeared the track 
energies, but not their directions, with a Gaussian error 



given by AE/E = 0.10/v 7 ^ + 0.01 (E in GeV) for lcp- 
tonic tracks and AE/E = 0.5/-\/i?-l-0.03 for all hadronic 
tracks. 

Displaced vertices at the LHC were identified requiring 
that the neutralino decays away from the primary vertex 
point, that is, outside an ellipsoid centered at the primary 
vertex 
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(5) 



3 We note that in the upper left dark region the stau is the LSP 
and in what follows we will not consider this region. 



where the z-axis is along the beam direction. To be con- 
servative we assumed the ellipsoid size to be five times 
the ATLAS expected precision in each direction for the 
semiconductor tracker [221 ] which are 5 xy = 20 /im and 
5 Z = 500 /im. To reconstruct the vertices we required 
that visible tracks coming from neutralino decays must 
have an intersection inside a sphere determined by the 
tracking detector resolution which we assumed to be 
10 /im [22j. Furthermore, we considered only the charged 
tracks inside the pseudo-rapidity region of |r/| < 2.5. 

Since the Higgs production in the LSP decay is char- 
acterized by the presence of two b-tagged jets we looked 
for events with at least one displaced vertex containing 
at least one jet tagged as a b-jet. In our analyses we 
considered a b-tagging efficiency up to 50%. 

In order to ensure that the detached vertex events 
are properly recorded we accepted only events that pass 
very simple trigger requirements. We further required the 
events to present an isolated electron (muon) with pt > 
20 (6) GeV, or the presence of a jet with px > 100 GeV, 
or missing transverse energy in excess of 100 GeV. 

For our analysis we have fixed tan/3 = 10, A = —100 
GeV and fi > 0. For this choice of parameters, the Higgs 
mass lies in the range 110 GeV < Mh < 120 GeV when 
we vary Mo and Mi/ 2 - Since we are only interested in 
detached jets coming from Higgs decays, we have further 
required that the jet-jet invariant mass is around the 
Higgs mass value. 

Within the SM framework displaced vertices originate 
from decays of long lived particles, like i?'s and r's, and 
consequently its visible decay products exhibit a rather 
small invariant mass. In contrast, in our BRpV model, 
the displaced vertices are associated to the LSP decay 
and will have in general a large invariant mass associated 
to them. Therefore, physical SM processes do not lead to 
sizeable backgrounds to the detached Higgs searches due 
to to large difference in the invariant mass of the visible 
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products. However, BRpV LSP decays into vZ are a 
potential source of background for the Higgs signal. 

As an illustration we show in Figure[2]the jet-jet invari- 
ant mass distribution of all displaced vertices exhibiting 
jets. As we can see, a cut on the invariant mass outside 
the range 100 GeV < Mi nv < 125 GeV eliminates a good 
fraction of supersymmctric backgrounds coming, for in- 
stance, from the neutralino decay into W and Z bosons as 
well as the three-body bbv channel. The physical back- 
ground can be further suppressed by requiring that at 
least one of the jets associated to the displaced vertex is 
tagged as a b jet. Moreover, these requirements ensure 
that SM backgrounds coming from the decay of long lived 
particles are also efficiently eliminated. There remain in- 
strumental backgrounds [23j which require a full detector 
simulation along the lines we have described above; this 
simulation is beyond the scope of the present work. 



In Figure [3] we show that almost all vertices containing 
b-jets come from neutralino decay via Higgs and that 
our invariant mass cut will eliminate the vZ background, 
while keeping a large fraction of the signal events. We 
checked that the events passing the LHC triggers and all 
the above cuts come from the signal events Xi ~ * v ^> with 
the physics background being negligible. 



In order to estimate the LHC reach for Higgs search 
coming from displaced vertex signal in BRpV-mSUGRA 
models we considered a few scenarios. In the optimistic 
analysis we assumed that there is no event coming from 
instrumental backgrounds or overlapping events and took 
the b-tagging efficiency to be 50%. In this case we re- 
quired that the signal must have more than 5 events since 
no background is expected and present our result in the 
Mi/2® Mo plane for integrated luminosities of 10 and 100 
fb _1 . We also considered three additional scenarios. In 
the first one we studied the impact of a lower b-tagging 
efficiency (30%) but we still assumed that the process is 
background free. In the second case we assumed that 
there are 5 backgrounds events originating from instru- 
mental errors and overlapping events and required a 5cr 
signal for a 50% b-tagging efficiency. Finally, in the last 
scenario we assumed the same background as in the last 
case, lowering however, the b-tagging efficiency down to 
30%. 
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Figure 2: Jet pair invariant mass distribution in GeV. The 
light blue (greyish) histogram stands for the background 
where the lightest neutralino decays via W and Z bosons and 
the other histogram stands for the channels where the lightest 
neutralino decays into bb pairs. 

IV. RESULTS 

In Figure [4] we depict the LHC discovery reach for the 
Higgs displaced vertex signal in our most optimistic sce- 
nario. The shaded (yellow) region at the bottom stands 
for points already excluded by direct LEP searches while 
the upper-left corner of the M 1 / 2 <8> Mo plane, the (red) 
shaded area, corresponds to the region where the stau 
is the LSP [l3j, and hence is not covered by the present 
analysis. The region around M 1 / 2 = 200 GeV has no 
signal due to the fact that the neutralino mass is smaller 
than the Higgs mass in it, therefore, being forbidden the 
two-body LSP decay into Higgs-neutrino pairs. 

From Fig. H one can see that the ATLAS and CMS 
experiments will be able to look for the signal up to 
Mi/2 ~ 700 (900) GeV for a LHC integrated luminosity 
of 10 (100) fb _1 . Notice that the LHC Higgs discovery 
potential is almost independent of Mq. For a fixed value 
of Mi/2 the LSP total production cross section decreases 
as Mo increases, however, the LSP branching ratio into 
Higgs-neutrino pairs increases with Mo, therefore, both 
effects tend to cancel and produce the observed behav- 
ior. Moreover, this figure also exhibits the average decay 
length of the neutralino, demonstrating that its decay 
takes place inside the vertex detector, ensuring a good 
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vertex reconstruction. 

We have also estimated the reach expected at LHCb 
for our Higgs search proposal. The hatched region in 
Fig. [4] indicates the LHCb reach for 10 fb" 1 . Due to 
the strong cut on the pseudo-rapidity required by this 
experiment the reach for 2 fb -1 is severely depleted and 
only a small region of the parameter space is covered, i.e., 
300 GeV < M 1/2 < 350 GcV and 200 GcV < M < 500 
GeV. 

Tagging b-jets emanating from a detached vertex is 
certainly a more intricate procedure, therefore, we also 
considered a lower b-tagging efficiency in our analyses. 
Figure [5] contains the reach of LHC for Higgs search us- 
ing a b-jet reconstruction efficiency of 30%, instead of 
50% used of Fig. [5J however, we still assumed that the 
search is background free. Comparing Figs. 0] and El one 
can see that the LHC reach in this second case is mildly 
affected by this change for an integrated luminosity of 10 
fb , while the changes are minute at higher integrated 
luminosities. 

A study of the instrumental backgrounds and the effect 
of overlapping events does require a full detector simu- 
lation, which is beyond the scope of this work. In order 
to assess the impact of existence of non-physical back- 
grounds we considered that these backgrounds give rise 
to 5 background events for both integrated luminosities 
used in our studies. In Figure [5] we present the 5cr LHC 
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Figure 3: Invariant mass distribution in GeV of the neutralino 
decaying into b-jet pairs separated into its several channels. 



Higgs discovery potential assuming a b-jet reconstruction 
efficiency of 50% and 5 background events. We can see 
from this figure that the existence of background events 
does lead to a substantial reduction of the LHC reach for 
Higgs in displaced vertices. 

In Fig. [7] we present the reach of LHC for Higgs search 
in a very pessimistic scenario that exhibits a lower b— 
jet reconstruction efficiency of 30%, as well as, the pres- 
ence of 5 background events. In this case we observe 
a more severe reduction of the LHC reach that is re- 
duced to Mi/2 = 600 GeV at most. This large depletion 
of the LHC search potential follows from the need of a 
large number of signal events to establish the signal given 
the fast decrease of the SUSY production cross section 
with increasing M~i/ 2 . In this sense, the 100 fb -1 case is 
more affected since the production cross section exhibits 
a steep decrease for M 1 / 2 > 700 GeV. 
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Figure 4: LHC reach for Higgs search in displaced vertices for 
the BRpV-mSUGRA model in the plane M1/2® Mo assuming 
tan/3 = 10, A = -100 GeV, and fj, > 0. The yellow stars 
(blue squares) represent the reach for an integrated luminosity 
of 10 (100) fb -1 while the hatched region corresponds to the 
reach of the LHCb experiment for an integrated luminosity 
of 10 fb" 1 . The (yellow) shaded region in the bottom stands 
for points excluded by direct LEP searches, while the (red) 
upper-left area represents a region where the stau is the LSP. 
Note that the black lines delimit different regimes of LSP 
decay length. 
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Figure 5: Same as Fig.[4]using a b-jet reconstruction efficiency 
of 30% with no background events. 



Figure 7: Same as Fig.|4]using a b-jet reconstruction efficiency 
of 30% and 5 background events. 
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Figure 6: Same as Fig. [4] using a b-jet reconstruction effi- 
ciency of 50% and assuming the existence 5 background events 
for both integrated iuminosities. 



troweak theory. We have given a quantitative analysis 
within the simplest minimal supergravity model with bi- 
linear breaking of R-parity, which accounts for the ob- 
served pattern of neutrino masses and mixings observed 
in current neutrino oscillation experiments. Similar vari- 
ant schemes can be envisaged where, for example, super- 
symmetry and/or electrowcak breaking is realized differ- 
ently. 

In an optimistic background free scenario the Higgs 
search in LSP decays can be carried out for LSP masses 
up to 300 (380) GeV for an integrated luminosity of 
10 (100) fb _1 . We showed that this result is robust 
against variations of the assumed b-tagging efficiencies. 
Notwithstanding, the results change drastically if instru- 
mental backgrounds are present. Assuming the existence 
of 5 background events reduces the LHC reach to LSP 
masses of 210 (250) GeV at the low (high) luminosity 
run. 



V. CONCLUSIONS 

In summary we have seen how the search for displaced 
vertices containing b-tagged jets at the LHC may not 
only provide evidence for supersymmetric particles but 
also lead to the discovery of the Higgs boson of the elec- 
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